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]Q '. Abstract 

The Kubo formula for the conductance of a mesoscopic system is analyzed semiclassically, 
yielding simple expressions for both weak localization and universal conductance 
fluctuations. In contrast to earlier work which dealt with times shorter than 0(log h^ 1 ), 



here longer times are taken to give the dominant contributions. For such long times, many 

X' 

distinct classical orbits may obey essentially the same initial and final conditions on positions 



and momenta, and the interference between pairs of such orbits is analyzed. Application 
to a chain of k classically ergodic scatterers connected in series gives — 1[1 — (fc + l)~ 2 ] for 
the weak localization correction to the zero-temperature dimensionless conductance, and 
^[1 — (fc + l) -4 ] for the variance of its fluctuations. These results interpolate between the 
well known ones of random scattering matrices for k — 1, and those of the one-dimensional 
diffusive wire for k — > oo. 
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I. INTRODUCTION 



Semiclassical ideas have been central to the understanding of transport effects in 
mesoscopic systems, specifically Weak Localization (WL) and Universal Conductance 
Fluctuations (UCF), from the outset However, they are usually used for handwaving 

arguments, with the actual calculation being done by resummation of perturbation theory 
expressions [ffl (Feynman diagrams). In recent years, fabrication of ballistic mesoscopic 
systems has become feasible, and it has been demonstrated that the chaotic or integrable 
nature of the classical dynamics in such systems is reflected in the quantum interference 
corrections to their transport properties f|. The perturbative approach with respect 
to the impurity potential is inapplicable to such systems, and indeed calculations using 
other theoretical tools such as Random Matrix Theory § (RMT) and the Nonlinear 
ex-Model |7|] have recently appeared. The recent progress in applications of the Semiclassical 
Approximation (SCA) to classically chaotic systems |||| gives rise to a hope that 
quantitative results for the mesoscopic transport effects could be obtained from it. This 
hope is realized below, but it turns out to be necessary to apply the SCA in a somewhat 
unorthodox manner. 

The reason for the difficulty in the semiclassical description of WL and UCF is obvious: 
both effects involve quantum interference corrections to the classical conductance, which 
are smaller in powers of H. Various authors have dealt with this difficulty in different ways. 
Some have followed the diagrammatic derivation quite closely, and used it for calibration 
of the magnitude of the effect |§. Others have limited their attention to effects such as 
coherent backscattering, where the quantum interference corrections appear in the leading 



order, and are as large as the classical result [ 1Q| . Still others have concentrated on the 



magnetic field dependence of UCF, and used RMT to calibrate the magnitude [pl] -[T3[|. All 
of these studies assume either diffusive or ergodic classical dynamics, and do not give general 
semiclassical expressions. 

Many of the recent analyses have used the Landauer formula or scattering approach, 
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rather than the bulk approach of the Kubo formula. The proof of the equivalence of these two 
formulae relies on the unitarity of the quantum mechanical evolution (current conservation). 
However, semiclassical evolution in classically chaotic systems is only approximately unitary, 
and therefore application of the SCA to the Kubo or Landauer formulae may lead to different 
results. The experience gained from the diagrammatic calculations WM shows that higher 



order corrections to the propagators (diagrams with Hikami boxes fllql ) are necessary for an 
evaluation of the conductance from the Landauer formula, but not when the Kubo formula 
is used. For this reason, we use the latter in the present work. 

A semiclassical analysis has already been developed for the Kubo formula in classically 
chaotic systems by Wilkinson |[16|| , with recent applications to transport in antidot arrays 



17fl . It was assumed there JT6[ that the relevant propagation times are shorter than the 
Ehrenfest time, £g ~ 0(log so that a coherent state or wave packet maintains its 

correspondence to a point in the classical phase space throughout its evolution. However, 
it has since become known that semiclassical expressions are applicable also to later times 
| |18| |. At such later times the region in classical phase space which evolves from an initially 
minimal-uncertainty wavepacket is stretched and folded by the chaotic dynamics into a very 
long and curved shape, which may intersect the region defined by some final wave packet 
several times. The Ehrenfest time thus marks the onset of interference between different 
classical orbits which correspond to these different intersections (called the mixing regime), 
and not the breakdown of the SCA. In the present work the focus is on times much longer 
than the Ehrenfest time, and compact expressions for the contribution of these interference 
terms to the conductance are developed. Interestingly, some of the terms (including the 
periodic orbit contributions studied by Wilkinson) retain the same form whether the orbits 
under consideration are short or long on the scale of whereas other terms are peculiar to 
times longer than £g and can not be found by studying the strict semiclassical limit h — > 0. 

The semiclassical expressions derived below exhibit three rather novel features: firstly, 
it will be assumed that the actions of the classical orbits involved are large, and that they 
thus contribute with essentially random phases and a statistical description is appropriate; 



secondly, interference effects between continuous families of classical paths appear explicitly; 
thirdly, the distribution of classical paths is described as having a continuous density in phase 
space, even when the initial position and momentum are given, i.e. it will be necessary 
to introduce a small amount of averaging over the initial (and final) conditions. Such a 
description is relevant not only for the average over an ensemble of similar mesoscopic 
systems which differ in their microscopic details (the disorder ensemble), but also for a 
single system at times longer than the Ehrenfest time. For example, in the expression for 
weak localization, orbits which start at some given position r and momentum p, and end 
after some time t at the same position r but at momentum — p will be needed. Rather than 
retaining only self retracing paths, which would be implied by strict classical mechanics, the 
results involve the density of orbits in phase space around the point (r, — p) (note that a 
self-retracing path is not related to a distinct path by time reversal, and is thus irrelevant 
for interference between pairs of time reversed paths). In incorporating these features we 
are following the approach of our previous work ||19|| , which pertains to density-of-states 
correlations in mesoscopic systems; however the need to account for contributions which do 
not strictly obey the initial and final classical conditions did not arise there. 

As in other applications of the semiclassical method, interactions are ignored (apart from 
a possibly self-consistent potential), and the electron fluid is considered as a degenerate 
thermal distribution of non-interacting particles moving in a mesoscopic sample. We have 
in mind a situation where the electrons' motion is mostly free, but occasionally they hit an 
impurity or the boundary of the system and are scattered. This corresponds (in the language 
of quantum chaos ||) to the motion of a 'billiard ball' which is scattered by some obstacles 
(see Fig. 1). In order to apply the SCA, it is assumed that all the classical dimensions, such as 
the mean free path or the radius of curvature of the obstacles, are much larger than the Fermi 
wavelength (i.e. h is small). The classical dynamics is assumed to be completely chaotic, 
with all orbits hyperbolically unstable (systems with integrable or intermediate dynamics 
require separate consideration). The number of dimensions or degrees of freedom, N, must 
thus be larger than 1. The exponentially large number of orbits in such systems justifies 
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our statistical approach. Possible complications associated with caustics etc. are ignored in 
the present work, and only the generic contributions are considered. In addition we assume 
that the spectrum is essentially continuous, i.e. that the single-particle level spacing A is 
much smaller than the other energy scales in the system, such as those determined by the 
broadening of the levels T due to the external leads, the temperature T or the frequency 
uj (even when those tend to zero). This also implies that there is no exponential Anderson 
localization in the system even in two or quasi-one dimensions, i.e. the system is smaller 
than the localization length, because AcEc where E c is the Thouless energy (E c — k/to 
with t E the time for an electron to traverse the system). 

The various contributions to the conductivity and its fluctuations will be expressed as 
integrals over the distribution function of classical orbits in the system. Given a Hamiltonian, 
H{r, p), an initial point in phase-space, (r, p), and a propagation time t, classical dynamics 
generates a final point which we denote by (r t , p t ). The unaveraged distribution of classical 
orbits is defined as a 5 function around this point: 

f(r',pi,t;r,p) = 8(r'-r t )8(pi-p t ). (1) 

It is convenient to use a distribution function f E limited to the energy hypersurface, by 
factoring out the energy conservation condition: 

f E (r', P ', t; r, P ) S(H(r', p')-H(r, p)) = f(r', p', t; r, p) . (2) 

Whereas these distribution functions describe the chaotic classical dynamics in intricate 
detail, they can be averaged over small ranges in initial and final conditions to give a 
smooth distribution which describes the evolution in a statistical sense. This averaging will 
be denoted by an overline: f E (r', p', t; r, p), with the range of averaging determined by %. 
For times t significantly longer than the Ehrenfest time t E , the distribution f E becomes 
independent of the specific details of the averaging procedure. For example, in an ergodic 
system f E becomes independent of initial and final conditions at long times. Apart from the 
distribution function f E , additional properties of the classical paths occasionally appear in 
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the expressions (e.g. the distribution of areas enclosed by the classical paths is relevant for 
the case of weak magnetic fields). 

Apart from the phase-space averaging denoted by the overline, angular brackets will 
be used in order to denote averaging over the semiclassical phases, which are assumed 
to be uncorrelated except for possible symmetries, e.g. (e lSa l h ^ = and ^e l ( Sa ~ Sf3 ^ h '^ = 
8 a ,p + b a T ,pi where S a and Sp are the actions of classical orbits a and (3, and a T denotes 
the orbit time-reversed to a. This assumes that time-reversal is the only symmetry in 
the system, and that the number of orbits is exponentially large so that the possibility of 
self-symmetric orbits, a = a T , may be ignored. This second averaging should be understood 
either as an averaging over the specific positions of the obstacles (the disorder ensemble 
average), or, for a single system, as an averaging over a range of possible values of the Fermi 
energy // (or a range of values of Planck's constant H) . The average quantum correction to 
the conductance will be denoted by (AG), and its variance by ((AG) 2 ) Var . The fact that 
these quantities can be written semiclassically as integrals over the distribution function 
Je, without reference to individual classical orbits and the plethora of their actions and 
amplitudes, is very useful in applications to classically chaotic systems fl9| , |20|j . 



The outline of the paper is as follows: in Sec. II the general SCA expression for the 
Kubo conductance is derived, and shown to reduce to the classical conductance when all 
interference effects are ignored. The appropriate choice of the electric field distribution 
is discussed. Weak Localization is analyzed in Sec. Ill, and Universal Conductance 



Fluctuations in Sec. IV. In these sections the general expressions in terms of /e(r', p', t; r, p) 
will be applied to diffusive and ergodic systems, for which simple expressions for /g are 
readily available. We concentrate on calculating the magnitude of the interference effects 
at zero temperature and magnetic field, in order not to repeat the considerations given 
in the previous semiclassical analyses p|,|^JT0HT2| . Extensions of the analysis, e.g. to finite 
temperatures, and applications to more complicated systems which consist of several cavities 
connected in series through ideal leads, will be considered in Sec. V, followed by a discussion 
in Sec. VI. A short description of this work is being published separately plfl . 
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II. THE KUBO FORMULA 



According to linear response theory, the real part of the conductivity tensor is given by 
(see e.g. appendix A of Ref. [Q]) 

where £ m and \m) are the single-particle eigenenergies and eigenstates respectively, p is the 
momentum operator, and Jfd{0 — V( ex P[M£ — /•*)] + 1) is the Fermi-Dirac distribution, 
with fi the chemical potential and (3 the inverse temperature (e and m are the electron 
charge and mass respectively, u is the frequency, and vol is the volume of the system; we are 
assuming a simple effective mass description of the electrons, but the derivation applies with 
minor modifications to systems with a non-spherical Fermi surface). This can be rewritten 
in terms of the single-particle propagator [e.g. 5(H — e{) = J2 m \ m )$(£,m — e i)( m \] as 



e h f 
(Jjk(uj) = — - — / dei de 
m 2 vol J 



huo + ei — €2 + iO huo 
J drdr' (7^7 (r'|«5(^- e2 )|r>^) (r\5(H - e 1 )\v') ) j , (4) 

where we have used pk = f dr |r) j-£r ( r | ; and ignored spin (a spin index summation should 
be understood with each spatial integration; we avoid the extra indices below by assuming 
spin degeneracy). 

For the purpose of a semiclassical analysis, it is convenient to define quantities which 
are bilinear in the quantum mechanical propagators, and to transform them into a form 
depending on one energy and one time variable, rather than two energy variables [19|. We 
are thus led to define the following "form-factor" : 

K(E,t;r,r') = J de e^Tr (j(r') 5(E + \e - H) }{v)S{E - \e - H)) , (5) 

where j(r) is the current density operator (for a given spin projection). Some of the basic 
properties of this form-factor are: (a) it is real; (b) it is symmetric under the interchange 
of all indices and the sign of the time variable, Kjk(E,t;r,r') = K^j[E, — t; r', r); and (c) 



if an external magnetic field is the only source of time-reversal symmetry breaking, then 
reversing the sign of this field (together with the sign of the spin indices) has the same effect 
on K_ as reversing the sign of t. As we will discuss only weak magnetic fields and use the 



SCA, we can rewrite the form-factor as 122 



K jk (E, t; r, r') * ^/ dee^ (j^jW(E + §e - £)|r>) (jj^l^ - Je - H)\*')) ■ 

(6) 

The conductivity is given in terms of this form-factor as follows: 

/r dv 1 
dE dt Ke F(E, t) I — p K{E, t; r, r) , (7) 

where 

F(E, t)= e«/» ' fFD(E-je)-f FD (E+je) 

Causality is reflected by the fact that replacing u by — uj is tantamount to taking the 
complex conjugate of F(E,t). In order to elucidate the structure of F(E,t), note that 
(fFD(E—^e) — fFD(E+^e)^j /e plays a role similar to 5(E — /j) (it is a box function of 
width e and height 1/e, smeared by the temperature which will be assumed small; as 
we will be interested in long times t, only small values of e are relevant). In fact, it 
is straightforward to show that / dE F(E,t) = e luJt 9(t) — ^S(t), where 9(t) is the step 
function, and the last term cancels with the diamagnetic term. The detailed form of F(E, t) 
becomes considerably simpler if only the time symmetric combination of its real part is 
needed, ffie (^F(E,t) +F(E,—t)j. Due to the properties of K_ mentioned above, indeed 
only this combination appears if one is interested in a symmetric integral over K_, such as 
the longitudinal conductance, or if the magnetic field vanishes (in the cases to be considered 
below both of these conditions are satisfied). The integration over e then becomes trivial 
due to a 5 (e — flu) factor, and in the limit of small frequencies we may write 

g(co - 0) ~ | o °° dt J dE (-f FD (E)) J ^ K(E, t- r, r') . (9) 



In the zero temperature limit to be considered below, the integration over E may be omitted, 
and its value is identified with the Fermi energy \i. Notice that the part of F(E,t) which 
is asymmetric in t and responsible for the Hall effect has been simplified in this analysis, 
which assumes that a is a Fermi surface property [^] . 

In the case of restricted geometries or non-diffusive systems, the conductance rather 
than the volume-averaged conductivity is the appropriate quantity to study. One defines 
the space-dependent conductivity ^(r, r') exactly as in Eqs. (0) and (||), except that the 
integrations over r and r' and the division by the volume are not performed. The current 
is given by j(r') = / dr cr(r, r') E(r), where E(r) is the electric field. The dissipative 
conductance of a sample of general geometry with two leads, at zero frequency, can be 
written as the dissipated power / dr'E(r') • j(r'), divided by the voltage V squared: 

G=^Jdrdr' E(r') g(r, r') E(r) . (10) 

Current conservation is expressed as V r • a(r, r') = V r ' • cx(r, r') = 0, in the absence of a 
magnetic field [V r • V r > • <x(r, r') = still holds even if a magnetic field is present]. This can 
be shown by integration by parts in Eq. ( [TO] ) (see, e.g., Ref. |TJ]]) to imply that the electric 
field E(r) need not be calculated self-consistently, and instead one can use any electric 
field distribution which gives the voltage V when integrated along any path connecting 
the two leads [the boundary conditions require that the components of E(r) and of cr(r, r') 
perpendicular to an insulating boundary vanish]. In fact, one may take different electric 
field distributions for the two E(r) factors in Eq. (|10"D — in the scattering approach of the 
Landauer formula the electric field is concentrated in the source lead for one factor, and 
in the drain for the other (another example is the case of multilead devices, for which the 
conductance Gy is a matrix, and the boundary conditions for the two electric field factors 
may be different). As already mentioned, the semiclassical expressions derived below for 
cr(r, r') do not necessarily obey current conservation, because higher order corrections in h 
are not included. Thus, the SCA expressions for the conductance do depend on the use of 
the actual electric field distribution in the sample, as discussed below. 
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The next step is to write down the semiclassical expression for K(E, t; r, r') of Eq. (^j). 
As a starting point we use the van Vleck formula: 



(r'|exp(-iih)|r) = £ A a e iSa/n 

a£{r,r';i} 



(11) 



where a is a discrete index which runs over all the classical paths that start at the point r 
and end after a time t at the point r'. The amplitude and action of the classical path a are 
given by 

1/2 



det 



1 dp c 



h dr' 



S a = J' (p dr - H dt) 



(12) 



where p a is the initial momentum of the path a and v a is the integer Maslov index (this 
and other factors of i may be ignored for the purposes of the present work). After Fourier 
transforming from time t to energy E, we have (see, e.g., 0) 



{r'\5(E - H)\r) 



*£{r,r';E} 



(13) 



where the index a counts classical paths of energy E, and the modified amplitude and action 
are given by 



A, 



det fl 9Po 



h dr' 



1/2 



1 dT n 



h dE 



1/2 



S a 



p dr . 



(14) 



The derivative appearing in this amplitude is taken at a constant duration of the orbit 
t = T a (it is also possible to re-express this amplitude in terms of an (N+l) x (JV+1) matrix 
of derivatives taken at constant energy E, but this will not be helpful here; N denotes 
the number of dimensions). Again, the Maslov index z> Q may be ignored, since only the 
magnitude of A a will be needed below. It is necessary to note the derivatives of the action: 

dS a m dS a dS n 



dE 



dr 



dr' 



Pa , 



(15) 



where T a is the duration of the classical path a, p a is the initial momentum mentioned 
above, and p' a is the final momentum. 
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Substituting Eq. (p~3|) in the expression for the form-factor, Eq. (H), and using the fact 
that only contributions from small values of e will be important (because t is integrated over 
a large range) to develop the action around the mean energy E, gives 

K(E,t;r,r') ~ / dee"**/* ^ A a A} p' a p^ e^+l^-^+I^A 

m a,/3e{r,r';S} 



4 /» E ^ P« P/> c*^-^ 5 (t - ) . (16) 

"' a,/3€{r,r';E} 



The 5 function over time should be understood to have a width determined by the higher 
order corrections in e. Inserting this result in Eqs. @ and ( |T0D will give the semiclassical 
approximation for the conductivity and the conductance. 

Before discussing quantum corrections to the conductivity, we observe how the classical 
results can be regained from this expression. To this end, all interference terms in Eq. ([IB]) 
are ignored, and only the 'diagonal' part of the double sum, a = (3, is retained: 

K D (E,t;r,r') = —^h £ \A a \ 2 p' aPa 5 (t - T a ) . (17) 

771 «€{r,r';B} 



In order to proceed, the amplitudes Eq. (|14D should be substituted here. As similar 
expressions will be used below, we note here the general form of a sum of this type: 

h £ \A a \ 2 5{t-T a ){. . .) a = \A a \ 2 5(E a -E)(...) a 

aG{r, r';.E} a£{r,r';i} 

= jdro^o 5 ^ H(ropo) _ E ^ 5{rQ _ r) 5{rt _ rl) (...) {ropo) 

= I dpE dp ' E ^ E< ^'' p ' t] r ' ^ ^ ' '' )(r ' p) ' ^ 

Again, the phase space point (r t , p t ) is that which evolves from the initial point (r , Po) by 



1 dT a 
h dE 



in \A r 



following the classical dynamics for a time t. In the first equality the factor 
was used to turn from a fixed energy representation to fixed time. In the second equality 
the sum over a was rewritten as an integral over the initial coordinate and momentum, 
using the factor \A a \ 2 = det (jj^jjr) (the integration over the initial position is trivial, and 
in fact superfluous at this stage). The result is an integral over the phase space energy 
hypersurface, where the properties of the individual paths which were denoted by the dots 
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are now identified by the initial coordinate and momentum of each path. Notice that all 
the determinants of derivatives which appeared in the amplitudes have been replaced by 
integrations over 5 functions p0[| , in such a way that allows for the introduction of the 
classical distribution function f E (r', p', t; r, p) in the last line, [integration over the energy 
surface is denoted by / dp E . . . = J dp 5 \ H(r, p) — E^j . . .]. 

Applying this trick to the diagonal approximation of the form factor, Eq. fll~7l), gives the 
classical contribution to the conductivity of Eq. (|): 
p 2 

cl e 

(at zero frequency and temperature). This may be rewritten as 



1 f°° f 

T^J dt J drdp E dr'dp' E pp f E (r ,p ,t;r, p) (19) 



p2 m 

v d = —v \ dt (pp') t (20) 
— m z Jo 



where v = h N\ ol J dr dp 5 \ji — H(r, p)j is the density of states (implicitly including the spin 
summation), and the last factor is a momentum correlator: 

Jdrdp E dr'dp' E . . . / B (r', p', t; r, p) 

( ' ' - }t ~ T^dV E (21) 

(it is the classical counterpart of the Fermi-surface correlator (p(0)p(t)) defined in ]2|). 
Eq. (|20|) is the classical contribution to the Kubo conductivity. Notice that the density of 
states contains a factor of h~ N , so that the quantum corrections to the conductance, of 
the order of e 2 /h, are small corrections to it (higher powers of H). In the present work, 
the semiclassical limit is considered with the Fermi momentum and the mobility (or other 
characterization of the scattering potential) taken as classical parameters, so that the density 
of electrons and the conductance become trivially H dependent. 

For diffusive motion, given some initial value of the momentum p, the average of p' at a 
time t shortly thereafter is equal to p multiplied by exp(— t/r), where r is the momentum 
relaxation (or transport) mean free time. Integrating over the directions of the initial 
momentum, the classical conductivity is found to be diagonal, a% = 5j t kCr D , and 

a D ^) = —^ e J^ dt cos(ut) exp(-t/r) v ^ (22) 
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where the frequency dependence has been restored (pp denotes the Fermi momentum). Since 

2 

the density of electrons n is equal to ^f^jl, this evaluates to 

aDt u \ = I^gfe T __ . (23) 

m 1 — itu 

which is just the Drude conductivity. 

The simple description of the momentum correlations used above is inadequate for the 
calculation of the space-dependent conductivity. In fact, it is known from diagrammatic 



theory that when all interference terms are ignored, it may be written (for u = 0) as [14 



of fc (r, r') ~ ^(0) [^(v-r') - V,V^(r, r')] . (24) 

Here 8(r— r') is a smeared 8 function of range equal to the mean free path I, which represents 
the short range part of the conductivity (analogous to the Chambers formula), and is due 
to paths which have not scattered at all. The scaled 'diffuson' d(r, r') represents the long 
range contributions of paths which have scattered at least once, and obeys the equation 
— V 2 d(r, r') = 8(r — r') with vanishing boundary conditions at the conducting leads and 
vanishing normal derivative at insulating boundaries. This ensures that cr°(r, r') conserves 
current. 

It is possible to render the long range part unimportant by using the classical electric 
field distribution in Eq. (|T0|) (for the standard rectangular geometry this is just a constant 
field). Indeed, if the condition VE = is valid, then it follows from the boundary conditions 
on E and d(r,r') that / drdr' Vj V k d(r, r') E(r) E(r') vanishes by integration by parts. In 
this case one may use a simple approximation to the space-dependent conductivity, keeping 
only the short range part [the first term in Eq. fl24|)l, without compromising the accuracy 
of the classical conductance G cl . It is emphasized here that this choice of the electric 
field represents the actual electric field in the sample, if it is interpreted as the gradient of 
the electro-chemical potential, rather than the externally applied perturbation which can 
be arbitrary. In other words, if charge neutrality is assumed (the chemical potential can 
not vary, and the self-consistent electric field is just that which will not cause any charge 
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perturbations), then the long range part of the conductivity can not contribute, because it 
represents the currents due to the gradient of the induced charge perturbations. 

This idea can be generalized to non-diffusive systems, such as a chaotic cavity (see 
Fig. 2). We will assume that the classical dynamics in such a cavity is not only ballistic, but 
also ergodic, so that the probability for an electron at any point r inside the cavity to leave 
through (or to have come from) the left (right) lead is proportional to its width Wl (Wr)- 
The self-consistent electrostatic potential is thus a constant within the cavity, and is equal 
to {WlVl + WrVr) I '(Wl+Wr), where Vl and Vr denote the potentials in the corresponding 
reservoirs. All of the potential drop occurs in the leads (or at the boundaries between the 
leads and the cavity or the reservoirs). Note that the velocity correlator of Eq. (^) is 
multiplied by the electric field in Eq. (|T0|) and integrated over time, in such a way that the 
contribution of a certain path [determined by its initial conditions (r, p)] to the conductivity 
has a simple interpretation: it is proportional to the potential difference between the initial 
point r and the reservoir which an electron with the given initial conditions will eventually 
reach. If the electron passes through an ergodic cavity on its way, and if the electrostatic 
potential is chosen self-consistently, then it is no longer necessary to integrate along the 
remainder of the path — the potential in that cavity is already the averaged potential 
of the reservoirs, weighted by the probability that the electron would leave through the 
corresponding lead. The long-range part of the electron paths, i.e. following the electron 
all the way to the reservoirs, thus becomes unimportant (in the diffusive case any elastic 
scattering event which randomizes the electron's direction of propagation plays the role of 
an ergodic cavity, in that it ends the short-range part of the propagation). 

According to this discussion, the precise form of the electric field in the leads does not 
matter. For speficity, we take the classical electric field E(r) to be constant over regions of 
size a in each lead. Combining Eq. ( |T0D with Eq. ([19]), the classical conductance may be 
written as 

° d = ^JWyl I drd PEdr'dp' E (p-E(r)) (p' • E(r')) f E {v', p', f; r, p) . (25) 
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The contribution of the short range part to the conductance thus becomes (the integration 
over f E (r',p',t) is trivial): 

G ^lE^/^^P)4co^), (26) 

where Vp is the Fermi velocity, AVi denotes the voltage drop over the corresponding lead, 6 is 
the angle between p and the direction of the lead, r(r, p) is the time that an electron starting 
at (r, p) spends in the region a of the electric field, and j i dr denotes integration over the 
electric field region in the lead i. The integration over the directions in / dp E , together with a 
factor of Vp \ cos(#)|, may be replaced by an integration over the transverse momentum / dp± 
and a summation over the two possible directions along the lead. When summed over these 
two directions, the free time r(r,p) gives ajvpj\ cos(#)|. The integration over the position 
along the lead gives a further factor of a, which thus cancels out as it should. The remaining 
integral gives the number of transverse channels in the lead = h~^ N ~^ J { dr± dp±, which is 
proportional to the width of the lead (in two dimensions ^ = 2WiPp/h, in three dimensions 
Wi is an area and gi = nWipp / h 2 ) . The final result is (ignoring spin) 



e 2 

rid _ e 
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9R \ 9l \ 

9l+ [ ; 9r 



6 - 9R9L . (27) 

h 9R + 9L 



s 9r + 9l) \9r + 9l 

This result corresponds to adding the resistances of the two ideal leads classically in series. 
Obviously, in order to reach this result with any other choice of the electric field factors, one 
would have to evaluate also the contributions of the long-range parts of f E (for example, 
the transmission which enters the Landauer formula is due to paths that start in one lead, 
scatter inside the cavity, and then leave through the other lead). 

The issue of an appropriate choice of the electric field becomes much more important 
for the quantum corrections to the conductivity, to be considered below. The reason is 
that it is quite hard to find a current conserving approximation for cr(r, r') which includes 
these quantum corrections, even in the diffusive case which is treated well by perturbation 
theory (see Ref. [J3J). When the self-consistent electric field configuration is used, current 
conservation and the long range part of the conductivity become less important. For 
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example, if quantum interference gives an enhanced probability to find an electron at some 
(r, p) at some time t, it becomes unnecessary to follow the propagation of that electron to 
possibly correlated momenta p' at later times . Loosely speaking, the missing electrons 
represented by the non-strictly-vanishing value of V • j(r) may be thought of as being 
re-injected into the system with a random direction, so that the self-consistent potential 
at r automatically takes care of their contribution. Unfortunately, the validity of this 
approach can only be strictly proven if one can write down an expression for the higher 
order corrections to er (r, r') , and show that they do not contribute when integrated with 
the electric field factors. While this is readily done in the diagrammatic analysis, it is not 
easily generalized to other, non-diffusive systems. In the following sections we proceed by 
analogy with the diffusive case, and calculate the quantum corrections to the conductivity of 
a general classically chaotic system, using the self-consistent electric field configuration and 
the leading order quantum corrections to a(r, r'). We return to this issue in the final section 
and show that this approach is indeed justifiable, at least for an ideally ergodic cavity. 



III. WEAK LOCALIZATION 

In this section the SCA is used to calculate the average of the quantum correction to the 
conductivity, i.e. the weak localization correction P,p6fl. We concentrate on the long-range 
part of the conductivity, i.e. on times t larger than ~ r — the short-range part does 
not have a weak localization correction. As advertised, the actions of the classical orbits in 
Eq. ( |16|) for the form factor K_(E, t; r, r') will be assumed random and uncorrelated, except 
if the two orbits are related by a symmetry. After averaging, only two types of contributions 
remain: the classical contribution (3 = a, and that of interference between time reversed 
orbits (3 = a T (it is assumed that time reversal symmetry is the only symmetry in the 
system). For any orbit a 6 {r, r';E}, we have a T G {r',r;E}, with a and a T sharing 
the same values of action, amplitude and duration, but p a r = —p' a and p' aT = — p a . The 
possibility of having a strict equality f3 = a T arises only in the case that r' = r, giving rise 
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to a factor of 2 enhancement of K^E, t; r, r) relative to its classical value. As the coordinates 
r and r' are integrated over, it is necessary to find how this enhancement is reduced when 
r' deviates from r. Therefore, we include in the weak localization term all pairs of orbits 
for which (3 ~ a T , in the sense that (3 and a smoothly deform into a pair of time reversed 
orbits when r' approaches r. Orbits which are self-symmetric are excluded, because their 
contribution is already accounted for in the classical term. 

The expression for the weak localization correction to the form-factor thus reads 

(AK(E,t;r,v')) ~ — 2 h £ \A a \ 2 6{t-T a ) p'^-p'J e^ § ^ n , (28) 

m Q,/3e{r,r';E} 

f3~a T yta 

where it is assumed that r' is near r, so that the only important r' dependence is in the rapidly 
varying phase factor. A direct evaluation of the r and r' integrations over the form-factor 
in the H — ► limit gives vanishing results. In fact, the stationary phase conditions would 
imply Pa = P/3 and p^ = p^, which in classical mechanics can only hold for self-symmetric 
orbits, a = (3. In principle, one could try to evaluate higher order corrections to the 
non-stationary phase integrals which arise, in the h — > limit. However, in practice (cf. 
Ref. ||) h is not extremely small, and the number of possible orbits a in a chaotic system can 
be exponentially large. In the mixing regime (t > ts) many of these orbits have such small 
momentum differences p^, — p' aT so as to make the phase practically stationary throughout 
the integration region. The spatial integration region is limited by the size of the system (in 
practice it may be smaller because the paths a and f3 may cease to exist due to caustics or 
shadowing). Thus, the contribution of orbits with momentum differences smaller than h/l± 
is just proportional to the size of the integration region, l±. Rather than following the exact 
distribution of possible values of l± for different orbits, in the following we describe this 
result effectively by a 5 function over the momentum difference. Note that this represents a 
non-standard application of the SCA, because it is assumed that many orbits can fit into the 
width of the 5 function, which is proportional to h. We return to this point in the discussion 
of Sec. VI A. 

The weak localization correction to the form factor, Eq. (|28|) , is again in the form 
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described by Eq. (0), which allows us to re-express it in terms of the classical distribution 
of orbits. Furthermore, Eq. (|T5D may be used to develop the actions around the point r' = r: 



S a (r,r';E) ~ S a (r, r; E) + (r'-r)p^ ; 

S p {r, r'; E) = S a (r>, r; E) ~ S a (r, r; E) - (r'-r)p a (29) 

(inclusion of higher order terms in this expansion turns out to be unnecessary), giving: 

(AK(E, t; r, r')) * ~ / rf P£ dp^ / s (r, p', f ; r, p) p' p' c ^-rW)/» . ( 30 ) 

In order to perform the integrations using the stationary phase approximation, it is necessary 
to treat the pre-exponential factor in the integrand as slowly varying. To this end, we 



replace /^(r', p', t; r, p) with f E (r',p',t;r,p), where the overline denotes averaging of the 
initial and final positions and momenta over small ranges. As long as the range of averaging 
is much smaller than a Fermi wavelength in position, and much smaller than a typical value 
of h/l± in momentum, there is no way that this replacement can affect the result of the 
integrals in Eq. ( |50D and Eq. (|S[). After the replacement, the pre-exponential factor is in 
fact smoothly varying, if the time t is indeed significantly longer than the Ehrenfest time 
ts [it is not necessary to explicitly remove the contributions of self-symmetric orbits from 



/e(r, p', t; r, p) — such orbits are exponentially rare in the mixing regime]. 

The next step uses the stationary phase approximation, specifically / dx dpf(x,p)e lxp ^ h ~ 
hf(0, 0) for small % and a smooth f(x,p), in order to perform the integrals over the angular 
variables of p' and the transverse components of r'. The stationary phase conditions identify 
p' with — p, and the components of r' perpendicular to p with those of r. The longitudinal 
component of p' is not integrated over because of the limitation to the Fermi surface, which 
gives rise to a factor of 1/vp. The integration over the longitudinal component of r', parallel 
to p, can not be done in the stationary phase approximation. However, as occurs also in the 



case of the spectral form-factor [T9|], this integration is trivial — the integrand is constant - 
and the result is just equal to the effective length of the integration region, which we denote 
by l(r, p). This integration region is not limited by a small H, and may extend over relatively 



long distances along the direction of the classical path: it represents constructive interference 
between a continuous family of classical orbits, labeled by the longitudinal component of r'. 
As a result, one has 

/ > e 2 1 f°° f Z(r p) 

(Act) = -- -/ dt / drdp E pp/ B (r, -p,t;r,p) — - — , (31) 

\ _/ m z h vol jo J vf 

which is a general semiclassical expression for the weak localization correction to the 
conductivity [the precise meaning of l(r,p) will be discussed further below]. Note that 
the integrations over r and r' always lead to an identification of p a with p^, and of p^ 
with p'p (cf. ||22||). In the present case of the weak localization contribution, we also have 
an identification of p with — p' (see Fig. 1), leading to the negative sign of the complete 
expression. 

For diffusive behavior with isotropic scattering, the meaning of /(r, p) is identified with 
the free path (see Fig. 1) — when r' deviates from r further than the next or the previous 
scattering event, the momentum factors in Eq. ( |T6D become essentially random (as opposed 
to the approximation used in Eq. (P5[)], leading on the average to a vanishing contribution. 
In N = 3 dimensions, and for times t > r, one may describe the distribution of classical 
orbits by 

Mr',p',t;r,p) = W(r, r'; , (32) 

where W(r, r'; t) dr' is the probability for a diffusing particle which started at r to be 
within dr' of r' at time t (in N = 2 dimensions 2np replaces 47rp 2 in the last factor). The 
factorization of the distribution into separate spatial and momentum space dependencies, 
with Je independent of the momentum direction, implies that the factors of pp may be 
replaced by dj^pp/N. The averaging over the integration segment l(r, p) gives 2vpT — the 
factor of 2 is due to Z(r, p) being defined as the sum of the free paths in the backward and the 
forward directions [another way to justify this factor is to recall that a classical path with 
a long integration segment l(r, p) for the integration over the longitudinal component of r', 
will also have the same integration region for that component of r, which means that in the 
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averaging l(r, p) is weighted by its own length]. The resulting correction to the conductivity 
is again a diagonal tensor: 

2,sD t T 4> r dv 

(Aa jk ) = -e 2 6 j:k — J dtj— W(t, r; t) , (33) 

where D = v 2 f t/N is the diffusion constant, and an s = 2 spin degeneracy factor has been 
explicitly restored. This result coincides with Eq. (3.8) of Chakravarty and Schmid |§, 
and thus their quantum mechanical derivation of the numerical prefactor in this equation 
(their Appendix D, which assumes a white-noise disorder potential) may be replaced by a 
semiclassical one. 

The time integration appearing here is limited from above by the dephasing time 
70, which is due to interactions of the electron with other particles not included in the 
single-particle Hamiltonian (the assumption of a continuous spectrum implies that the times 
involved are shorter than h/A). It is also limited from below, by the mean free time r. The 
probability density W satisfies the diffusion equation: 

dW{T dt r ' ] t} - DWlWiv, r'; t) = 6(t)5(r - r') . (34) 

For short times, W(r,r';t) = (47rDt)~ N / 2 exp (— |r — r'\ 2 /ADtj. For times of the order of 
the diffusion time through the sample to = L 2 /D, it is necessary to expand W in the 
eigenfunctions $ n of the diffusion operator, W(r,r';t) = J2 n ^n( r )^n( r ') ex P( — t/T n ). In 
this representation the integrals over W(r, r; t) in Eq. ([33|) become trivial, and for — > oo, 

e 2 2sD 

w = -^ r - (35) 



n 



For example, for a quasi-one-dimensional wire extending from x = to x = L (with 
finite cross sectional area, i.e. many transverse modes) the longest lasting eigenmodes are 
<& n (x) oc sin(7ma;/L) (where n = 1,2,...) with decay times of T~ x = 7T 2 n 2 D/L 2 . The 
boundary conditions are essential in determining this — they are closed in the transverse 
directions but open in the direction along the wire. This reflects the fact that trajectories 
hitting the latter boundaries will continue through the hypothesized ideal leads into the 
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reservoirs, and will not return to any point r' in the sample (/ W dt is essentially the 'diffuson' 
of the diagrammatic technique, which was mentioned earlier). Defining the dimensionless 
conductance per spin direction g, and using the fact that S n ( 7rn ) _2 = 1/6, gives (for the 
j = k = 1 component) 

= ; (A9> = 4- (36) 

which is a well known result for the diffusive wire geometry. 



In order to apply the semiclassical expression, Eq. (plf) , to a more general system, 
one must specify the meaning of the free path factor l(r, p), or in other words one must 
perform the integration over the longitudinal component of r' with care. Note that in a 
two-dimensional system the longitudinal direction may be defined as the locus of points for 
which the action difference S a (r, r'; E) — Sp(r, r'; E) vanishes — there is then no reason to 
stop the integration at the next scattering event (which in itself may be ill-defined for a 
smoothly varying potential). Following the discussion of the previous section for the classical 
conductance of such systems, the value of the integral is identified as the potential difference 
Ay(r,p), which is defined as the difference between the potential at the point eventually 
reached by an electron at (r, p) and the point from which it emerged. Strictly speaking, 
AV(r, p) is equal to either ±V or 0, because both the original and the eventual points are 
in one of the reservoirs. However, the quantities in Eq. ([H]) are averaged, and so we are 
led to define Z(r, p) = AV"(r, p)pf/ (E(r) ■ p) [the averaging in AV(r, p) is over the same 
range as in J'e ] ■ The integration over the longitudinal component of r' is thus effectively 
limited by the "Ehrenfest length" vptE- 

With this definition, the application of Eq. (|31|) [cf. also Eq. (|T0D1 to the chaotic cavity 



of Fig. 2, involves AV(r, p) = ±AV{, due to the fact that the escape time from the cavity 
is assumed to be much longer than £g. With the specific choice of the electric field within 
the leads as before, |E(r)| = AVi/a, this gives 

(AG C ) = —r £ tj^ / dt drdp E v F cos 2 (£) f E (r, -p, t; r, p) /(r, p) , (37) 
n „_£ v a Jo J % 
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where the factors of a cancel due to the longitudinal component of the dr integration, and 
the relation /(r, p)| cos(#)| = a. The time integral of the distribution function f& is found by 
requiring that the total number of electrons escaping from the cavity, which can be written 
as an integral over any crossection in the leads, J2%=r,l fi dp' ± dt /^(r', p', t; r, p), is 
equal to unity (the escape velocity vp\ cos(#) | is used as before to transform from f t dp B . . . to 
Jj dp± . . .). The time integral / Je dt is independent of the detailed initial and final positions 
and momenta due to ergodicity in the cavity, and is thus equal to h l ~ N /(gL + gn) (as long 
as r and r' are in the leads, p is in the inward direction, and p' is in the outward direction). 
The remaining integral gives just h N ~ 1 gi, so that 



e 2 



(AG C ) = -- 



9R Y 9l ( 9l V gn 



jR + gLj gR + gL \gR + gh gR + gh 



e g L g R f v 
h [gR + giY 



In the case of equal leads, g^ = g^, one obtains (Age) = —\ (with G = g^r), in agreement 
with RMT results 0. Applications to additional systems will be considered in Sec. V. 



IV. UNIVERSAL CONDUCTANCE FLUCTUATIONS 

We now turn to the off-diagonal terms (3 ^ a,a T , which do not contribute to g_ on 
the average, and calculate the typical magnitude or variance of the fluctuations. A close 
inspection of the perturbative derivation (Refs. PJ2"7|] and references therein) shows that 
there are in fact three different types of contributions. A sketch of these three types of 
diagrams and the corresponding classical paths is given in Fig. 3, and the semiclassical 
analysis will be detailed in this section. 

The contribution of all (3 ^ a, ot T paths to the conductivity of Eqs. @ and ( fL6| ) is 
denoted here by Aa_ ND . At first sight it would seem that each term, defined by a specific 
choice of a and /3, has an uncorrelated phase (S a — Sp) /h, and thus its absolute magnitude 
squared gives an independent contribution to ( |Aa| 2 ) . This is indeed true for terms with 

\ — / Var 

T a ,T/3 ~ tEi and forms the first type of contribution [Fig. 3(a) and subsection A below]. 
If one of T a , Tg is negative [note that according to Eq. ([TED their sum must be positive], 



22 



then the integrations over r and r' in Eq. (|9|) effectively 'join' the paths a and f3 into a 
single periodic orbit of duration \T a \ + \Tp\ (the path with negative duration —\T\ may be 
considered as starting at r' and ending at r after a time |T|). Since a single periodic orbit can 
be bisected into two segments in many different ways, all contributing with the same phase, 
this type of contribution can show significant interference between different a, (3 pairs. It is 
useful to first add up all the contributions for each periodic orbit, which will be labeled 7, 
and then consider the contribution to the variance from each such term. This is done below, 
and forms the second type of contribution [Fig. 3(b) and subsection B]. As a result of the 
momentum factors, and the fact that the integral of the momentum along a periodic orbit 
must vanish, these terms do not contribute in the simple cases considered here. 

The third type of contribution arises when T a or T@ is positive, but smaller than the 
Ehrenfest time. For definiteness, take < Tg < in this case the path a forms a periodic 
orbit, which returns to its starting coordinate and momentum after a time T a — Tp and 
then continues along the same direction as (3 for a time Tp. Thus, the action along this 
last segment cancels in the expression (S a — Sp)/h, and just as in the second case described 
above, this contribution is close to a periodic orbit. These periodic orbit contributions form 
the third and last type of contribution to ^|A^t| 2 ^ v [see Fig. 3(c) and subsection C]. 

In the following subsections, the semiclassical expressions for each of these three types of 
contributions are derived. For simplicity, we consider only the case for which the temperature 
and the frequency approach zero, and for which time reversal symmetry holds. 

A. Contributions with T a ,Tp ~ 

The different contributions to ^Acrj^ Aa^ 2 J?*^ will be denoted by F 1 , F 2 and F 3 , with 
the ji/ci, indices suppressed in most of the equations (the complex conjugate of the real 
quantity Aerj^ is taken for convenience in notation; it amounts only to exchanging the 
a and (3 indices). According to Eqs. (|9|) and (|T6|), the expression for lAer^l 2 involves an 
integration over four different spatial coordinates. The first type of contribution, F 1 (or at 
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least that part of it that will not vanish after averaging) comes from regions where both r 
coordinates and both r' coordinates are close to each other, and both copies of a and of (3 
coincide. Due to time reversal symmetry, there is also a similar contribution from the case in 
which these coordinates are interchanged. As will become evident, these two cases contribute 
terms with indices 5j 1 j 2 5fci,fe 2 an d Sj 1 ,k 2 ^k 1 ,j 2 respectively, which are otherwise identical. The 
first of these, denoted by F la gives 



2 i \ 2 
e/ 1 



oo roo 







^ a = 1^2 Vol ' ' ^ ' ^ ' dr+dV '+ I dV - dr> - 



h 2 E KH^| 2 P«P.PLp^(ti-^)^(t 2 -^)e^. (39) 

The summation here includes only orbits with positive times T a and Tp, which scatter at 
least once, and are not related to each other by symmetry (other contributions are included 
in the other terms). The notation r + = |(ri+r2), r = ri — r2 etc. is used, and only the 
contribution for which ot\ = a.2 and (3\ = $2 is retained. The deviation of the orbits from 
{r + ,r^,yu} leads to corrections to the actions, which to first order are given by 

AS ~ -r_(p Q - p p ) + v'_(p' a - p' p ) . (40) 

Consider first the integrations over the time variables. Necessarily t\ = t2, and this 
time variable will be denoted by t + . It will be convenient to add a fictitious integration 
over / dt- 5 (t- — {T a — Tp)^j. The 5 functions over time may then be rewritten in the form 
S(t a — T a )6(tb — Tp), where t a = t + + t-/2 and tf, = t + — t-/2. After this is done, it is 
possible to transform the sums over a, (3 G {r + ,r' + ; fi} into phase space integrations using 
again Eq. fllSP- This gives 



,2 ^ \ 2 



/■'" = [ ^2^jj J t dt a J t dt b J dr + dr' + j dr_dr'_-^ J dp aE dp' aE J dp bE dp' bE 



t E 

/fi(r'+, P' a , ta] r+, Pa) f E (r' + , p' b , t b ; r+, p b ) p' a p b p' a p b e iAS/h , (41) 

where the quantities relating to the possible a orbits are denoted by a subscript a, and those 
of the (3 orbits by b. It is assumed here that the relevant contributions come from orbits 
longer than the Ehrenfest time, so that t a and t b are bigger than t E - This allows us to 
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replace the two factors of by their smooth averages and to neglect the contributions 
of a — (3 to these averaged distributions. The next step is to use the phase factor in order 
to perform the integrations over the relative coordinates and momenta, giving 



F la ~ (±-J—\ [°° dt a ^ dt b [ drdp E [ dr'dp 
\ m z h vol I Jo Jo J J 



E 



f E (r>, p', t a , r, p) f E (r>, p', t b ; r, p) ^ p) f ' P ' } p'pp'p, (42) 

where again the integrations over the longitudinal directions give factors of the 'free path' 
/(r, p), just as in the evaluation of the weak localization term. 

For diffusive behavior, Eq. ( |52"D can again be used, and the momentum direction 
integrations performed. The free paths Z(r , Po) are replaced as before by factors of 2v e t, 
giving 

x 2 /2e 2 L 2 \ 2 

r n^r^rJ dr dT ' W(r > r ' ta) W(r > r ' h) (43) 

(again the spin degeneracy factor s = 2 has been restored). Using the decomposition 
of W(r, r'; t) in terms of orthonormal eigenfunctions, the spatial integrals give simply 
J2 n exp(— t a /T n ) exp(—tb/T n ). Together with the term F lb , this yields 

F l = (^^) 2 (S jld2 5 klM + 5 jlM S kl>j2 ) E (^) 2 . (44) 

This expression can be seen to coincide with the first part of Eq. (46) of Altshuler and 
Shklovskii p7 |. For the example of the quasi-lD wire, with T n D/L 2 = l/ir 2 n 2 , this 



contribution to the conductance fluctuations gives (|Ag| 2 ) Var = ^ X^i n ~ 4 — i n terms 
of the dimensionless conductance g. 

A further application of Eq. (^) is to the chaotic cavity of Fig. 2. This gives (including 
the factor of 2 due to time reversal symmetry, but ignoring spin) 
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/ 2\ 2 AV 2 AV 2 r°° 



t E J tE J i J j 



f E (r', p', t a ; r, p) f E (r', p', t b ; r, p) /(r, p) l(r , p ) cos (0) cos (0 ) v F 

/ 2 \ 2 n 2 « 2 

2 (-) E 



^/ iJ=L,R (9r + 9l) 4 (9r + 9l) 2 
2 fgV , glg2 \ A . (45) 



Here we have used our previous result for the time integral of f E , and all the factors of a 
have canceled as before. The notation i denotes the lead opposite to the lead i. This is just 
twice the square of the weak localization result, and in the case of symmetric leads ga = gi 
reduces to the well known result ((Ag c ) 2 ) Var = ^ (see Ref. [§]). Notice that in this case 
all of the conductance fluctuations originate from the first type of contribution, Fig. 3(a), 
because there are no periodic orbits which traverse the region in which the classical electric 
field does not vanish. 



B. Contributions with T a T p < 

Consider next the contributions which are concentrated around periodic orbits. For 
definiteness, assume that Tp is negative. The description A*^e~ iSl3 ^ h of the path (3 G {r, r'; E} 
may then be replaced by the identical term ib^e* 5 ' 3 '^, associated with the path f3' G {r', r; E} 
where (3 corresponds to retracing f3' backward in time (Tp = —T@ > 0). One may rewrite 
Eq. (H) as 

K(E, t- r, r') ~ 4 h £ A a A p , p' a p'p e '<*«+V>/» 5 (t - ^> ) , (46) 

0'e{r',r;E} 

which is completely equivalent, but more convenient if the times T a and Tp are positive. 
When the spatial integrations over the form factor are performed, it is seen that indeed the 
stationary phase points occur when ppi = p' a and p'p, = p a . This means that the path (3' 
must continue the path a, and vice versa — together they form a periodic orbit. 
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The periodic orbits of energy E may be enumerated by the discrete index 7 G 
{r = r', p = p'; E}. Each periodic orbit is in fact a continuous family of periodic classical 
trajectories, which differ from each other by the choice of the initial position along the orbit. 
The contributions of all the different pairs of paths a and (3 which fall along the periodic 
orbit 7 must now be found. Note that if the integral over all positive times t is taken as 
indicated in Eq. (H), the last S function may be replaced by a restriction to T a > Tpi. The 
next step is to undo the Fourier transforms which led to the energy representation in terms 
of A a and Ap>, and to return to a representation in terms of A a and Apr. 

g2 I , ,00 ntl 



roc r r roc rt\ 

/ dt / drdr' K(E,t;r,r') - dr dr' / dt x / dt 2 e iEtl/h e 

Jo J — m z h J Jo Jo 



iEt 2 /h 



A a Ap,p' a p'p l e i ^ +s ^ n . (47) 



{r,r';t^ } 
/3'e{r', r; t 2 } 



The relation ~ (instead of ~) is used to indicate that only the contribution with Tp < is 
included here. The contribution of a periodic orbit 7 may now be calculated. The spatial 
integrations here resemble the convolution formula for the propagator at time ti + t 2 in terms 
of the propagators at times t\ and t 2 , of which the trace is then taken. Thus, apart from the 
integration over the time difference t\ — t 2 and the appearance of the momentum factors, 
the stationary phase integrations can be performed just as is normally done for the periodic 
orbits involved in the Gutzwiller trace formula PJj|,pO||. This gives 

dt J dr dr' E_(E, t; r, r') ~ 

— 2 J2 APe i§ ^ n ( Tl dt / T7/2 ^ 2 p 7 (t +ti)p 7 (to) , (48) 

™\ €{r=r f; =p , E} Jo Jo 

where p 7 (t) denotes the momentum along the periodic orbit 7, at a point parameterized by 
a time variable t, and t 1 +t 2 = T T The periodic orbit amplitude, 



1^1 



h 



-1/2 

det[M 7 -J] , (49) 



where M 7 is the monodromy matrix describing the stability of the orbit 7, is the same as 
the Gutzwiller amplitude apart from a time factor. As in the case of the non-periodic orbits 
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of the propagator, the squares of these amplitudes may be written as integrals over the 



distribution of classical orbits [[LJJ. 

\AP\ 2 5(t-T,) = j^- J^drdpE f E (r,p,t',T,p) , (50) 

where the integration is restricted to the region in phase-space surrounding the periodic orbit 
7 (the time variable too is restricted — the right hand side contains additional contributions 
at times t = mT 7 with any integer m). 

Eq. fl48| ) implies in fact that these periodic orbit contributions vanish in all the cases 
considered in the present work. One may extend the time integral of Eq. ([|) to negative 
times, because one is interested only in the contribution which is symmetric with respect to 
reversal of the magnetic field ||23|| . Another way to make this point is to recall that the time 
reversed orbit 7 T will also contribute, with p 7 t(£) = — p 7 (— t). Thus both time arguments 
appearing in the momenta p 7 (t) can be taken to vary over the whole periodic orbit. However, 
the factors of J 1 dtp y (t) must vanish for a periodic orbit, in order for it to return to the 
starting point. Such a cancelation has also been observed on the diagrammatic level PS 



where the condition T a ,Tpr > corresponds to the condition that a pair of propagators 
has either both propagators advanced or both retarded. This kind of contribution will not 
vanish only if the frequency u is non-zero, or if the integral J E • dr along the periodic orbit 
is non-zero (this can happen if the electric field is driven by a time dependent flux, and 7 
surrounds the flux). These exceptions will not be discussed further here. 



C. Contributions with T a or Ta ~ 



We now turn to the last type of contribution, which involves cases when either \T a \ or 
\Tp\ is smaller than the Ehrenfest time t^, so that the stationary phase conditions force the 
two orbits to overlap over this period. In this case too the integrations over the end points r 
and r' lead to the appearance of periodic orbits. In fact, most of the analysis of the previous 
subsection still holds, but with Tp < 0: the particle is performing periodic motion along a 



28 



path of period T 7 = T a — Tp (assuming for the moment that Tp < t E ), and after revolving 
around for a time T a it retraces its path for a short period Tp, and reaches its starting point. 
The corresponding contribution is 

/*oo r 

J dt J dr dr' K(E, t; r, r') ~ 

— 2 A?e^^ h [ Tj dt f dt 2 p 7 (t -t 2 )p 7 (to) • (51) 

7G{r=r',p=p';£/} a 

Note that the only difference between this and Eq. (|48| ) is in the limits of the t 2 integration. 
The absolute square of each such term is an independent contribution to the variance of the 
conductance, which thus includes a sum over all periodic orbits 7. 

Assuming that the relevant periodic orbits are longer than the Ehrenfest time, T 7 > t E , 
we identify the dt2 integration in Eq. (|51"D as an integration over the momentum correlator 
along the orbit 7. A periodic-orbit dependent diffusion constant may be defined, such 
that the time integrations in Eq. fl51"D give simply m 2 T^D_^. In order to express the results 
in terms of the phase-space distribution function fg, this diffusion coefficient is relabled as 
D_(r,p), which is identical to for all points (r, p) on the periodic orbit 7. Using this 



notation and the expression for the amplitudes, Eq. (p0[ ), gives 

/*oo r 

dtt J dr dp E f E (r, p, t; r, p) J2(r, p) I?(r, p) . (52) 



2 \ 2 



A factor of 2 arises due to contributions with T a < t E rather than Tp < t E [these give the 
complex conjugate of Eq. (|5"ID1, and a further factor of 2 allows for time reversal symmetry, 
i.e. contributions of pairs of orbits with S~, = SLt. 

In the case of diffusive motion (with s = 2), and assuming t ^> r ~ t E , the diffusion 
constants D_ can be approximated by their average D5j t k, giving: 



-(3. 

where (see Ref. Hl9| ) 

I 



fijiM $hM I d osc (n) I , (53) 



\d osc (fi)\ 2 = s 2 — dttjdr W(r, r; t) . (54) 
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In these expressions one sees most explicitly the observation made by Altshuler and 



Shklovskii |27j that the corresponding contribution to the UCF are associated with the 
fluctuations in the density of states (d osc ). The integrations may be done explicitly, giving 
again a simple sum over the extinction times of the modes: 

F3 - 4s2 (^) 2 ^.^?(¥) 2 - < 55 > 

This is identical with Eq. ([BP apart from the spatial indices, again in agreement with 
Eq. (46) of Altshuler and Shklovskii. It enhances the fluctuations of the conductance by 
50%, which thus totals for our example of the quasi-lD wire to ((Ag) 2 ) Var = ^| = ^, once 
again a well known result. 

For the calculation of the conductance of chaotic cavities, Eq. (|10D, one must take into 
account the fact that the electric field factors are position dependent, and may vary over 
the region covered by a periodic orbit 7. The periodic-orbit-dependent diffusion constants 
should be taken to reflect this: 



D, = f Tl dt f dt 2 (p 7 (t„-*2) • E(r 7 (t -* 2 ))) (p 7 (t ) • E(r 7 (t ))) . (56) 



The limit — tE of the dt% integration is to be understood in the same manner as the integration 
over the longitudinal component of the r' integrations above it is useful to define a 



free time r(r, p) = A + V(r,p)m/ (E(r) ■ pj, in analogy with the free path /(r,p) defined 
above. The potential difference A + V(r, p) is the difference between the reservoir eventually 
reached by electrons starting at (r, p) and the potential at r (the integration is only over the 
forward direction, not the backward one [|29|). Again, although the propagation times to the 
reservoirs may be long, the averaging effectively limits the length of the path contributing 



to r(r, p) by the Ehrenfest length. Strictly speaking, the averaging in A + V(r, p) should be 
taken over the region in phase-space corresponding to a minimal-uncertainty wavepacket 
defined by the "width" of the periodic orbit 7, i.e. by the Monodromy matrix M 7 . Thus, the 
fact that the classical path for electrons starting precisely on a periodic orbit never reaches 
the reservoirs is irrelevant. Our assumption that the Ehrenfest time is much smaller than 
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the typical propagation times in the system (such as T 7 for the relevant orbits) means that 
A + V(r, p) becomes a smooth function already with a significantly smaller averaging region, 
and therefore the precise form of the averaging is unimportant. 
With this definition, Eq. ( |56D becomes 

D, = F dt (p 7 (t ) • E(r 7 (t ))) 2 r(r,p) . (57) 

Notice that despite the averaging involved in defining r(r, p), the resulting D(r, p) function 
can in principle have large fluctuations, so that it becomes important to use the actual 
fluctuating /^(r, p, t; r, p) in Eq. (|52[) , rather than the smooth averaged one. A possible 
way to avoid this is to use a weighted average for D(r, p), i.e. define an averaging such that 
D 2 (r, p, t) = f E (r, p, t; r, p)D 2 (r, p) / f E (r, p, t; r, p). 

As already mentioned, the periodic orbits do not contribute to conductance fluctuations 
in the case of the chaotic cavity of Fig. 2 because the electric field vanishes in the region 
in which periodic motion can occur, so that D(r,p) = 0. A simple example of a class of 
chaotic systems for which the periodic orbits of Eq. (|52"D do contribute will be considered in 
the next section. Note that the relationship between the contribution of periodic orbits to 
conductance fluctuations and their contribution to the fluctuations in the density of states 
is essentially modified by the presence of the electric field factors, and no longer follows 
Eq. (H). 

V. APPLICATIONS TO SPECIFIC SYSTEMS 

In the previous sections, the semiclassical formulae for weak localization and universal 
conductance fluctuations were applied to two simple systems, the diffusive wire and the 
ergodic chaotic cavity. In the case that the leads are of a constant width, both of the quantum 
interference effects are stronger for the diffusive system than for the chaotic scatterer (— | 
vs. —\ for weak localization, ^ vs. ^ for the variance of the conductance). Furthermore, 
introducing assymetry in the chaotic cavity by taking gi ^ 9r decreases the interference. In 
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the present section systems which are intermediate between these two cases are discussed: 
in the first subsection a system of two cavities connected in series is considered, allowing 
for any combination of widths of the different leads; in the second subsection a string of k 
chaotic cavities in series is considered, with all leads equal in width. In the course of the 
treatment of these systems, results which are valid more generally, for networks of ergodic 
cavities connected by ideal leads, will be given. 

Before embarking on the detailed treatment of these special systems, several possible 
extensions will be mentioned. Consider first the effects of symmetry braking by weak 
magnetic fields or by spin-flip or spin-orbit scattering. This requires knowledge of the area 



distribution associated with the paths contributing to /^(r', p', t; ,r,p), or the distribution 
of 2 x 2 spin scattering matrices along them, and has been considered in Ref. @ for 
weak localization. Because of the close parralellism between the present approach and the 
diagrammatic theory, it is not surprising that the results for diffusive systems are reproduced. 
Specifically, in the case of complete symmetry braking the weak localization correction either 
vanishes, if time reversal invariance is broken, or is multiplied by a factor of — ^ for the case 
of symplectic symmetry (strong spin-orbit scattering). The variance of the conductance is 
reduced by a factor of 2 or 4, respectively. These results hold whether the system is diffusive, 
ergodic, or simply chaotic. 

Temperature and frequency dependences may likewise be treated, and again the 
diagrammatic results will be reproduced for the diffusive case. Observation of Eq. fll~6|) shows 
that the frequency couples to the sum of the periods of the two interfering orbits T a + Tg, 
whereas temperature (averaging of E in a range around /i) limits the contributions from 
orbits with large period differences \T a — Tp\. This gives different behaviors for the four types 
of interference effects considered (weak localization, and the three different contributions 
to conductance fluctuations). The temperature and magnetic field dependences of UCF, 
for ergodic cavities without time-reversal symmetry, were very recently studied by Efetov 
|7]] using the nonlinear cr-model. As noted there, the results are in agreement with 
previous semiclassical analyses HTUHl^l, except for the amplitude of the effect which has 
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been reproduced semiclassically only in the present work. The most striking result of this 
reference is that temperature smearing, as opposed to dephasing, does not change the typical 
area which enters into the magnetic field dependence, although it does change the form of 
that dependence somewhat. This is of direct relevance to the experimental work ||30|| , as 
it may allow the measurement of the dephasing rate from the magnetic field dependence, 
by inferring this typical area. It may be intuitively explained by noting that as the finite 
temperature limits only the difference between T a and T 13 , the areas they encircle may be 
arbitrary, and the area difference will be of the same order as the typical area which enters 
in the zero temperature case. Dephasing, on the other hand, limits the sum of T a and T", 
leaving only the contributions with typically smaller and differences. 



A. Two ergodic cavities in series 

The system considered here consists of two ergodic cavities connected in series through a 
lead with qm conducting modes, with gi and gn denoting the number of modes in the right 
and left leads as before. In this case ergodicity is not achieved, and there are two different 



regions of initial conditions which must be considered in order to evaluate /^(r', p', t; , r, p). 
For electrons originating in the left cavity, the probability of leaving the system through 
the left lead is the sum of a geometric series: Pll = — — ( 1 — 



9L+9M V 9L+9M 9R+9M 

— 9L< - 9R+9K ^ whereas the probability of leaving the system through the right lead is 

9m9r~\~9r9l~\~9l9m ' o j & & 

P LR = — 2m as— fl 2M sm_ V 1 = 9m9r _ The probabilities for an 

^'^ gL+9M 9R+9M V 9L+9M 9R+gM J 9 M 9 R~\~ 9 R9 L~\~ 9 L9 AI ^ 

electron in the right cavity (or situated in one of the leads but headed towards the 
right cavity) is given by similar expressions for Prl and Pr^r. The classical electrostatic 
potential in the left cavity is therefore 9l ^ 9r+9m ^ Vl+9m9rVr an d that in the right cavity 

^ J gM9R+gR9L+9L9M ' ° J 

9r(9l+9m)Vr+9m9lVl 

9m9r 



. The potential drops in the left, middle and right leads are 

9M9R+9R9L+9L9M & 



AVr = V 



9m9r + 9r9l + 9l9m ' 



AV M = V- ' //, " // ' 



9m9r + gRgh + 9l9m 
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V 



9l9m 



(58) 



9m9r + 9r9l + 9l9m 



The classical conductance for this system may be written as 



( 



) 



2 



e 2 9l9m9r 

h 9m9r + 9r9l + 9l9m 



H 9i 



V 



(59) 



i=L,M,R 



This is the result of adding classically in series the resistances h/(e 2 gi) of the three leads [cf. 
Eq. ( p7|) for the case of a single cavity] . 

The weak localization correction for this system is given as a sum of four terms, 
corresponding to different (r, p) integration regions in Eq. (|31"D: electrons in the left 
lead moving towards the left cavity, in the middle lead moving towards either the left 
or the right cavity, and in the right lead moving into the right cavity. It is convenient 
to generalize the probabilities P^j mentioned above, so as to allow indices which describe 
electrons in the middle lead, moving either to the right % = M_^, or to the left % = 
Pij = J dt Jj dr' ± dp' ± fE(r', p', t; r, p), where (r, p) can be any phase space point in the 
'directed lead' i, and the time integration excludes very short times for which no chaotic 
scattering has occurred. The sixteen 'probabilities' Pij can be obtained from the four 
probabilities P{j with i,j = L, R by noting that i = is equivalent to i = L, i = 
is equivalent to i = R, and similar equivalences can be obtained for the final condition, j, 
up to factors of 9m/ 9l and 9ul 9R-, respectively (the directions of the arrows for i = L,R is 
obvious and omitted in the notation). In analogy with Eq. (|3~8"D, one obtains 



where the index i T denotes motion time reversed to that denoted by i. In the case of leads 
of equal width g L = g M = g R , this gives (Ag 2 c) = —8/27, which is intermediate between 
the —1/3 result for diffusive systems, and the —1/4 result for completely random scattering. 

The F 1 and F 3 types of contributions to the fluctuations of the conductance may also 
be evaluated in a quite straightforward manner. The first gives, as in Eq. (|45|): 




9l9m9r{9l + 9m)(9m + 9r)(9r + 9l) 
h {9l9m + 9m9r + 9r9l) 3 



(60) 
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' 4 / i,j=L,M^,M->,R y 3 



i -J 



' 2\ 2 



2 ( "H 9l9m9r{9l + 9M)(gR + 9m) 



x (gL + gMKgfi + gMKgl + g|) + ^9l9r9m ^ 
{9m9r + 9r9l + 9l9m) % 



Again, in the case gi = 9m = 9r we find ((Ag 2 c) 2 ) 1 = ^fjj, which is intermediate between 
the 2/16 and 4/45 results of chaotic and diffusive systems respectively. 

In order to find the contribution of periodic orbits from Eq. (0), consider a periodic orbit 
which traverses the middle lead n 7 times, which must be even. The integral of Eq. ( |5TD 
gives _D 7 = (n 7 /T 7 ) (AV^-/2) [the one half comes from r(r,p), which is on the average 
a/2\ cos(9)\vf for r in the electric field region]. Eq. (|52"D then gives 

<( AG ^ 2 )L = 4 (x) 2 j _ M ^ M ^r^/^dp ± Mr,p,t;r,p)^^, (62) 

where we have replaced the integration over the whole length of the periodic orbits, 
Jq°° dt J drdpE fE{r,P,t;r,p) (n(r, p, t) 2 /4tj by an equivalent term which integrates only 



over the cross-section of the lead, J °° dt J2i li dr± dp± /e(i", p, t; r, p) (n(r, p, t)/A^j . The 
notation n(r, p, t) identifies the value of n 7 corresponding to the different periodic orbits, as 
was done with D(r, p) above [ |3"T| . 

Due to the appearance of the n(r, p, t) factors, this result cannot be written directly 
in terms of the 'probabilities' Pjj, and requires instead the following modification: the 
geometric sum in Pm^,m<_ and Pm^,m^ is weighted by a factor of m/2 for the mth term in 
the sum (n 7 = 2m), and thus gives \z/{l — z) 2 where z = g 2 M / '((?l + 5 , m)(5 i m + gR.) ■ Taking 
into account the contributions from both and M_+, one has 

«^> 2 >i - 4 iff L M+ :z+ g J w^ ■ < 63 » 

which for g L = g M = g R (with z = 1/4) gives ((Ag 2 c) 2 ) 3 = y^- This too is intermediate 
between the zero result for the single chaotic cavity, and the result of the diffusive wire. 
Adding the results of Eqs. ( |5T] ) and ( |6lf ) gives for the total conductance fluctuations 
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^(AG , 2c) 2 ) Var = 2 (~^) 9 2 l9m9r(9l + 9m)(9r + 9m)(9l + 9r) 



2\ 2 



x (g£ + flMXflfl + 0m)(0l + ga) - ^9l9r9m ^ 
(9 m 9r + 9r9l + 9l9mT 

(\ 5 
I J , which is again intermediate 

between the and results. 

16 15 

The quantum interference effects obtained in Eqs. ( |60| ) and fl6~j|) can be considered in 
the two limiting cases of a very wide or very narrow middle lead. In the case Qm 9l,9r, 
the middle lead connects the two cavities very efficiently, so that a particle in one of them 
will explore the phase-space of both cavities ergodically before having a chance to leave 
through one of the external leads. Thus, the results of a single cavity, Eqs. (38) and (45), 
are obtained in this limit. This can be understood already from the electric field distribution, 
Eq. fl58f ), which vanishes in the middle lead in this case. In the opposite limit gu <C 9l,9r 
the middle lead acts as a weak link (the use of semiclassics assumes g>j 3> 1, so it is still much 
larger than a quantum point contact), and the quantum interference effects are suppressed, 
giving (AG 2C ) * -(eV/OfltfGfc 1 + g R l ) and ((AG 2C ) 2 ) Var * 2(e 2 /hfgUg^ + g^f. As 
the expressions of Eqs. (|60|) and (|6~3]) are symmetric with respect to any permutation of the 
indices L, M, R, it does not matter whether the middle lead or any one of the other leads 
is taken as very wide or very narrow. However, the ratio of the periodic orbit contribution, 
Eq. (|6T|), to that of interference between pairs of different paths, Eq. (|6"3D, is not invariant 
under such permutations. In fact, this ratio is maximal when gu *C gi = 9r, and in that 
case it is equal to 1, which is larger than the ratio of 1/2 familiar from diffusive systems (it 
is for a single cavity). The quantum interference effects turn out to be largest when the 
three leads are of equal width, which is the case to be studied and generalized in the next 
subsection. 
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B. Chains of ergodic cavities 



It is clear that by adding more and more ergodic cavities in series, the situation of the 
diffusive one-dimensional wire may be approached. Consider the case of k cavities connected 
by leads of equal width, g>j = const. Even though the leads may be identical to each other, 
the chaotic cavities are taken to be different in order to avoid a periodic situation in which 
Bloch states would emerge (i.e. we assume the absence of translation symmetry). In this 
case the electric field configuration is trivial, with AVi/V = l/(k + l) in all leads. In order 
to make use of the semiclassical formulae of the type appearing in Eqs. (0), ( |6T| ) and (|63|). 
one needs to study and generalize the classical probabilities Pjj. 

It is convenient to define a classical dynamic probability pi tTn (t) equal to the probability 
that an electron will be found in the mth cavity, given that it was initially in the Zth cavity 
and that it has since traversed through a lead t times (l,m = 1, . . . , k and t = 0, 1, . . .). 
Notice that the 'time' variable t is discretized, with no reference to the actual time the 
electron may spend in the leads and in the cavities on its way. Quite generally, Pi j = 
jj L ~ J2tZoPi,m{t), where I is the cavity that the 'directed lead' i is flowing into, m is the cavity 
out of which j is flowing, and g m is the sum of the conductances of the two leads connected 
to the cavity m. The simple electric field configuration mentioned above is directly related 
to these 'probabilities'. 

A dynamic difference equation, similar to the diffusion equation, may be written for 
Pi,m(t): 

Pl, m (t+ 1) = \ [pi, m -l{t) +Pl,m+l{t)] , (65) 

with the initial condition pi >m (0) = Si }Jn and the boundary conditions pi,o(t) = pi,k+i(t) = 0. 
The classical evolution can be decomposed into i — 1, . . . , k eigenvalues a, and eigenfunctions 
Pi t i, such that Pi jfn (t) = Y$=i Pi,i<Xi Pi,m- F° r the simple system considered here explicit 
expressions are available: ctj = cos ^j7r, and @ij = J sin ^rpYvr. Just as has happened for 
the diffusive case, the eigenfunctions turn out to be unimportant, and only the eigenvalues 
— 1 < a,i < 1 will play a role. 
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Generalizing the results of Eqs. (0), (0) and fl63|) to the present case, and rewriting 
them in terms of the dynamical probabilities, gives for the weak localization correction 

{AGkc) = -JJkTW^ pM = -T(T^?i^ = -U( 1 -(fcTiy) ' (66) 

where the explicit values of were used in the last equality. The result for the conductance 
fluctuations is 

= (tt) (Fny( 2 ?(i- a ,p + ?(i-' a .,)2) 

= (tt) p^ T jl|([(*+ 1 ) 2 - 1 ][ 2 (*+ 1 ) 2 +7] + [(*+l) 2 -l][(*+l) 2 -4]) 

-(ffM 1 -^?)- (67) 

where the two different contributions are kept separate until the last equality. Notice that 
the factor of | = ^ in the relationship between P it j and pi >m exactly compensates the fact 
that the number of terms in the i,j summations is twice larger than in the /, m summations 
(there are twice as many directed leads as cavities). 

Naturally, these results reproduce those described above for k — 1 and k = 2 cavities 
(k — describes an ideal wire), and approach the diffusive results for k — » oo. A similar 
crossover behavior has been calculated for a different kind of system which interpolates 
between the same ergodic and diffusive limits, using supersymmetry techniques |32|]. In that 
case too the idea was to connect k ideally chaotic cavities in series, but the connections were 
made directly in the Hamiltonian, by introducing matrix elements which mix between the 
quantum states of adjacent cavities (each cavity was ascribed a GOE Hamiltonian). The 
results were as here rational functions of k which reproduce the k = 1 and the k — > oo limits, 
but their form was much more complicated than Eqs. ( po] ) and (|B7]) . Even the asymptotic 
behavior for large k is not similar, and in fact the crossover for UCF described in Ref. |32| 
is slightly non-monotonous. It would be very interesting to compare our results also with 
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a continuous crossover between an ideal lead (k = here) and a diffusive wire, obtained 
when a disordered region of size L and mean-free-path I is introduced in the lead (with 
the continuous parameter L/l). In principle such an analysis can be carried out using the 
semiclassical methods developed here, but would depend on the details of the disordered 
wire (Je and E(r) depend on whether the scattering is isotropic or small-angle scattering). 



VI. DISCUSSION 

The main results of the present work are the semiclassical expressions for the mean and 
variance of the quantum interference corrections to the conductance, i.e. weak localization 
and universal conductance fluctuations. For classically chaotic systems, these quantities are 



expressed as integrals over the distribution of classical orbits, /e(r', p', t; r, p), involving also 
additional quantities which can be derived from this distribution: the self-consistent electric 
field E(r), and the effective free paths Z(r, p) and diffusion constants D(r, p). Knowledge 
of this distribution of classical orbits is readily available in the applications considered here 
(see Sec. V), due to the assumption of a diffusive system or a system consisting of several 
ergodic cavities connected through ideal leads. The important task of demonstrating the use 
of these expressions on a generic system lies beyond the scope of the present work. However, 
it is stressed that finding this classical distribution numerically for a given potential should 
be relatively easy, because only a statistical knowledge of the classical orbits is necessary, 
and there is no need to form a full database consisting of exponentially many orbits. 

It was originally thought U that the SCA could help to bridge the gap between the 
down-to-earth experimentalists and the abstract mathematical analysis of the theorists. 
Unfortunately, the recent developments in the theory have considerably widened this 
gap, with the introduction of diffusion equations in n-dimensional eigenvalue spaces, and 
supersymmetric techniques. It is hoped that the present method will contribute to reversing 



this trend P3fl, although it is acknowledged that the semiclassical approach has its own 



limitations (in particular, the case of a few, or partially open, channels is outside its 
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scope). The following two subsections give a detailed discussion of the novel features of 
the present work, and in that context an attempt is made to bridge a different gap — that 
between the theory of disordered systems and the theory of quantum chaos. Possibilities for 
cross-fertilization between these two fields, in both directions, are pointed out. 



A. Modified Semiclassical Approximation for the Mixing Regime 

It should be emphasized that the semiclassical analysis was applied here in an unorthodox 
manner. A strict stationary phase argument would require the initial momenta (and final 
momenta) of the two classical paths a and f3 appearing in the semiclassical expression for 
the conductivity, Eq. (|16D, to be strictly equal to each other. As these two paths start at the 
same position, they would have to be either identical to each other (the classical contribution 
to the conductivity), or to differ by completing a different number of revolutions around a 
strictly periodic orbit [IB] [as in F 3 of Eq. (|52|)1. Thus no weak localization corrections, 



and no contributions of the first kind [F 1 of Eq. (|39|)] to universal conductance fluctuations 
would be obtained. However, the value of h is never infinitesimally small, and thus the 
stationary phase argument is never completely strict. In fact, in chaotic systems interesting 
contributions to physical quantities often arise from the mixing regime, i.e. from orbits with 
propagation times between the Ehrenfest time and the typical escape time t esc . This 
regime disappears in the extreme h — > limit, and it is thus not surprising that its detailed 
treatment deviates from the standard SCA. 

Consider for example the classical orbit of Fig. 1, and its behavior when r' is not strictly 
equal to r. Labeling the impurities by digits, we refer to the interference term between the 
path a = r-l-2-3-4-2-l-r' and the path (3 = r-l-2-4-3-2-l-r', in which the order of scattering 
has been reversed. Clearly S a and Sp are guaranteed to be equal only if r = r'. It is also clear 
[from Eq. (|T5|)1 that when r' deviates from r, the difference of actions will be proportional 
to the difference in the final momenta of these two paths, which for the case depicted in the 
figure is about 1% of pp. This momentum difference, p^, — p^, is always perpendicular to p' a 
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and never vanishes, so there is no strictly stationary phase contribution to the integration 
over the perpendicular component of r'. Furthermore, the magnitude of this momentum 
difference is roughly independent of the perpendicular component of r', as it is determined 
by the first part of these paths, which is identical for a and (5 (in the present example the 
first two scatterers are identical). However, the range of this integration, which we denote 
by l±, is finite — roughly 10% of the size of the system, L, in the example — and thus 
the action difference will never grow beyond l±(p' a — p'p) (both these factors depend also 
on the parallel component of r', but their product is roughly constant). Now the question 
of whether this pair of orbits contributes or not depends on the value of H: if it is much 
smaller than 10~ 3 p F L there will be a fluctuating (never stationary) phase factor, and the 
contribution will be unimportant, while if h is much larger than 10~ 3 p F L the phase will be 
negligible throughout the whole integration range, and there will be a finite contribution. 

Quantitatively, such contributions may be described by a factor of h N ~ l 5±(p' a — p^), 
where the 5 function is over the perpendicular components of the momentum, and has a 
finite width ~ h/l± and height ~ In the analysis of Sees. Ill and IV the averaged 

distribution function f E was introduced already at an earlier stage and the integration over 
r' and p' was done in one step, so that this S± function never appeared explicitly. For 
long times there are exponentially many classical orbits from r to r', with initial and final 
momenta within h/l± of p and p' respectively, and thus the averaging in / E (r', p', i; r, p) 
gives a smooth distribution. The fact that f E is smooth makes it invariant to the details of 
the averaging procedure. Thus there is no need to keep track of the actual distribution of 
sizes of l± (at least as long as the main contributions come from long orbits). 

The crossover between 'short' and 'long' times for this purpose occurs at an Ehrenfest 
time t E , which depends on H and the typical Lyapunov exponents. A pair of paths starting 
very near to each other in phase space will have their momentum difference multiplied by 
a factor of order l/R after each collision with an obstacle, where / is the free path and R 
is the radius of curvature of the obstacle. Thus, if the initial momentum difference is h/l± 
and the momentum difference at t E is required to be of order p E , the resulting estimate for 

41 



the Ehrenfest time is 

\og(p F l ± /h) 

tE ~ r n^/7*r (68) 

(suitably averaged values of r, / and l± should be used here). Often this estimate can be 
taken to imply t# ~ r, especially in diffusive systems for which it is customary to take the 
limit of small scatterers R <C I (recall that diffusive motion is not a valid description for 
times t ~ r anyway). However, in the semiclassical limit H will become so small that we will 
have tE t, and in principle this should by taken into account. For instance, the magnetic 
field dependence is determined by the distribution of areas enclosed by diffusive paths of 
length t, which should be replaced by t — 2t E since the paths are so close to each other that 
they enclose a negligible area for a time t E near their beginning and end (this argument 
does not apply to the contribution of periodic orbits). Since the dominant contributions 
come from long times, such corrections are unimportant for diffusive systems, except for 
exceedingly small values of H. They are also unimportant in the chaotic systems studied 
here, because of the assumption of ergodicity which was used for the individual cavities 
(we have essentially assumed that both ts and the time taken to traverse each cavity are 
much smaller than the escape time from that cavity, which typifies the length of the shortest 
relevant orbits). It would be very interesting to study systems for which is not negligible 
compared to the typical propagation times. 

The physics associated with the Ehrenfest time can be clarified by considering a different 
geometry. It is often claimed in the context of disordered systems, on the basis of a 
perturbative analysis, that any disordered system can be described by a single parameter 
- the diffusion constant — on all length scales larger than some physical cutoff such as the 
transport mean free path or the grain size for a granular material. However, an additional 
length-scale related to the Ehrenfest time occurs naturally in the semiclassical analysis. 
Furthermore, this length-scale depends not only on the microscopic characteristics of the 
potential, but also on h, and diverges for 7i — > 0. As an example, consider a thick slab 
of transparent disordered material. If it is illuminated by a well-collimated beam of light, 
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a speckle pattern may be observed in the transmitted light. Based on the perturbative 
analysis, one may expect that the condition on the thickness of the slab is given by the 
transport mean-free-path, but in fact it is clear that the relevant length is related to the 
Ehrenfest time (provided that the microscopic features of the disordered sample are 'soft', 
i.e. they must be larger than a wavelength and must not act as beam-splitters). For slabs 
of thickness intermediate between these two lengths, the beam will remain well collimated 
and no speckle pattern will be observed, even though the direction of propagation of the 
transmitted (or reflected) light will be random. 

It is somewhat surprising that the importance of the averaging in /^(r', p', t; r, p), or 
equivalently the width of the 5±(p' a — p' /3 ) functions, was not emphasized earlier. This may 
be due in part to the fact that much of the attention was devoted to physical quantities 
which involve periodic classical orbits, such as the density of states |8|f|. In fact, this 
kind of averaging does not arise naturally in the analysis of the contribution of periodic 
orbits to the conductance either (see Sec. IV C above), and even if it were to occur, it 
would have no dramatic effect on integrals involving /^(r, p, £; r, p). In contrast, in the case 
of weak localization the averaging is around the point /#(r, — p, £; r, p), and it introduces 
qualitatively new types of orbits (i.e. non-self-retracing orbits) into the calculation for long 
times. Likewise, the averaging for /e(r', p', t a \ r, p) /^(r', p', t&; r, p) has an important effect 
as it allows distinct classical orbits to overlap [see Fig. 3(a)]. 

B. Corrections to the Conductance From Leading Order Propagators 

It is remarkable that the use of the Kubo formula, with a self-consistent choice of the 
electric field, has enabled here the calculation of corrections to the classical conductance 
which are of higher order in h, without having to evaluate such high-order corrections to 
the propagators themselves. As explained already in Sec. II, the choice of the electric fields 
can be shown to be unimportant on the basis of unitarity, and it is somewhat surprising that 
the semiclassical results for g_(r, r') do not obey unitarity, order by order in %. However, it 
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turns out that the leading order results for £(r, r') may give sub-leading contributions to the 
derivatives V r ■ a(r, r'). In order to clarify these issues, we briefly reconsider the example of 
a cavity with ideally 'random' scattering, which can be described by the circular orthogonal 
ensemble (COE) of RMT ||. As described in the following paragraphs, this system can 
be observed to display all of the above surprising aspects (i.e. a non-unitary leading order 
approximation which still reproduces the quantum corrections to the conductivity correctly, 
if the Kubo formula is used rather than the Landauer formula), without invoking any 
semiclassical considerations. 

Take the n x n scattering matrix, S^-, which is associated with such a 'random' 
cavity, connected to leads with a total of gi, + gR = n conducting channels. The 
dimensionless conductance is given in terms of this scattering matrix by the Landauer 
formula gc = J2ftiYl]=g L +i \^i,j\ 2 - Taking S it j to be a random member of the COE, one 
may calculate various averages of its matrix elements for any n. For example, the fact 
that (\Sij\ 2 ) = (1 + Sij) I (n + 1) implies that (g c ) = giga/in + !)■ Expansion in 1/n 
gives (gc) — gLgain 1 — n~ 2 ), in agreement with the semiclassical results for an ergodic 
cavity, Eqs. ([27]) and (p8|). Evaluation of the variance of the conductance is slightly more 
complicated, although straightforward 

((9c)% av = ^((|5 M f) 2 ) Var + ^(n-2)((|^| 2 )(|^ fc | 2 )) Var 

+9L9R(9L9R-n + l) ((|5 M f)(|^| 2 )) Var 

= 9L9R (^3) " (^TIf) + 9l9r ^-^ (^3) " (nTiy) 

. / n + 2 1 \ 

+9L9R{g L gR-n+l) [ n{n+1){n+3) ~ ^y 2 ) 

* ^WrK • (69) 

Here the indices i, j, k, I which appear in the different type of covariance terms are all taken 
to be different from each other. The result is again in agreement with the semiclassical 
result of Eq. (f4"5[) for an ergodic cavity. Note that it is obtained due to a cancelation of 
the first two contributions with each other, and that knowledge of moments of Sij to the 
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second sub-leading order in 1/n turns out to be necessary. Obviously, if the leading order 
expressions for the moments of Sij were used in the Landauer formula, only the classical 
conductance could have been calculated correctly. 

Consider now the freedom of using the known property of unitarity of the scattering 
matrices, J2j l*Si,j| 2 = l for each row or column. This allows us to replace the contribution of 
each row of the scattering matrix to the Landauer formula, Y^j= 9L +i \Si,j\ 2 by t ne combination 
(1— a)(l — Z)f=i \Si,j\ 2 ) + OiYJj =9L+ i \Si,j\ 2 , with any value of a. As the self-consistent classical 
electric field for an ergodic cavity is concentrated in the leads (cf. Sec. II), it is possible to 
follow a line of derivation equivalent to that used above for the SCA simply by choosing an 
appropriate value of a, specifically a = gi/n. Repeating this replacement for each row and 
each column gives the "Kubo formula" for the conductance, 

9o = - £ BkEjlSuf , (70) 

with the "classical electric field" factors 

f g R /n iil<i<g L 
Ei={ (71) 
I -9L/n if 9l+1 <i<n 

(the relative minus sign is due to the field being directed into the cavity in one lead, and 

out of it in the other). Observe that the classical conductance of the cavity is given here 

by a "short-range part", which does not require knowledge of the "long-range" scattering 

matrix at all. It is stressed that Eq. (|70|) follows directly from the Landauer formula and the 

unitarity of Sij. In the present context it is considered to be more basic. Evaluation of the 

mean quantum correction to the conductance can now be performed using only the leading 

order result for the corresponding moment of the matrix elements, (\Si t j\ 2 ) ~ (l + 5i j)/n, 

despite the fact that in this approximation the scattering matrix is not unitary. Most 

terms cancel with each other because J27=i^i = 0> an d the remaining — Y^i=\E 2 /n gives 

the WL result. The results for UCF may likewise be obtained using only the leading order 

expression for the moments, {(\Sij\ 2 ) (\Sk,i\ 2 )) Yar — (8i t k8j,i + 8i,i8j,k)/n 2 . The analogy with 

the semiclassical calculation is thus complete. 
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It is expected that explicitly using the self-consistent electric field could be very helpful 
in other calculation schemes too. For example, it should enable the calculations of Ref. |7j 
to be done by expanding only to the second order in the diffuson-Cooperon expansion, and 
not to the sixth order as was found necessary there (this is based on counting the number 
of ladders in the different diagrams for the "long-range part" of UCF ) . 



In this system one can explicitly study the next order corrections to the moments of 
Si j, and observe how their contributions could in principle be as large as the WL and UCF 
terms being calculated, but vanish due to cancellations between the different terms. This 
cancelation can be demonstrated without actually calculating the higher order corrections, 
by arguing that YU Ei — an d using the fact that the COE is insensitive to the actual values 
of the % and j indices. As noted in the closing paragraph of Sec. II, a similar cancelation can 
be demonstrated for diffusive systems using the diagrammatic technique. More generally, 
one may argue on the basis of universality that additional contributions from higher order 
terms should not be expected, at least in the cases of systems involving ideally ergodic 
cavities (Sec. V). In order to verify this expectation, it would be very interesting to study 
explicitly the contribution of higher order corrections which are known for semiclassical 
propagators [33J], especially those which are known to involve Ti 1 ^ 2 or h 1 ^ 3 corrections, i.e. 



the effects of diffraction and caustics. 

In the field of disordered systems, higher order corrections in h are usually analyzed in 
terms of two distinct small parameters . The first gives the accuracy of the description of 



individual scattering events, and may be written as Xf/1 where is the Fermi wavelength 
and / is a microscopic length such as the mean free path (for a semiclassical analysis, an 
obvious microscopic length is the radius of curvature of the obstacles depicted in Fig. 1). 
It is often argued that inaccuracy with respect to this parameter can be tolerated, because 
the actual potential in a real system is unknown anyway, and instead the mobility or the 
effective scattering cross-section are measured directly. The second small parameter is 
essentially h/At esc ~ 1/g, where A is the mean single-particle level spacing, t esc is the 
length of time a typical electron spends in the sample before leaving through the leads, and 
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g is the dimensionless conductance of the sample. This parameter describes the extent of lack 
of unitarity in the semiclassical approach as discussed above. Recently, it has been suggested 



that unitarity could be built into the semiclassical approximation from the outset ||36|| . It 
would be very interesting to analyze this approach from the point of view of two distinct 
small parameters. Such an analysis could imply that the SCA becomes, when unitarity is 
enforced, analogous to the nonlinear cx-model — the latter describes a disordered system 
with a white-noise potential to all orders in 1/g, but only to leading order in 
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FIGURES 



Fig. 1: Sketch of a path which together with its time reverse contributes to weak 
localization. The initial and final points may be varied along the segment /. The source (S) 
and drain (D) regions are connected through ideal leads to particle reservoirs. 

Fig. 2: Schematic sketch of conductance through a ballistic cavity, which is considered to 
be completely chaotic and ergodic. The electric field is taken to be concentrated in regions 
of size a in the left and right leads (the results are independent of a). 

Fig. 3: Sketch of the three different types of interference contributions to the conductivity, 
together with some of the corresponding perturbative diagrams. The classical paths a and 
P start at r and end at r'. (a) a different from f3, (b) a and (3 lie on a periodic orbit 7, 
and one of them is of negative duration, (c) a and (5 lie on 7, and (5 is short (in this case a 
traverses the whole of 7 and then repeats the segment (5 a second time). The sketch assumes 
diffusive motion with t E ~ r. 
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